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Abstract 



We discuss the enumeration of Feynman diagrams at tree order for processes 
with external lines of different types. We show how this can be done by 
iterating algebraic Schwinger-Dyson equations. Asymptotic estimates for 
very many external lines are derived. Applications include QED, QCD and 
scalar QED, and the asymptotic estimates are numerically confronted with 
the exact results. 
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1 Introduction 



With the growing complexity of scattering amphtudes that are becoming 
amenable to calculation, especially with the availability of recursive algo- 
rithms like im, 0, the question of the number of contributing Feynman graphs 
becomes of interest in its own right. As will be seen, already for tree diagrams 
their enumeration is a nontrivial question. In |^ the number of graphs for 
theories with a single, self-interacting scalar field was studied, and a method 
derived to estimate the asymptotic number of tree graphs for the 1 — > n am- 
plitude for large n was described, improving on earlier estimates 0, |^. The 
precise number of graphs, but not its asymptotic form, for the case of QCD 
was discussed also in 0, albeit in a manner that does not lead to straight- 
forward numerical results. It is the aim of the present paper to improve on 
this situation. In section 2, we briefly review the single-field case. In section 
3, we discuss the case of 'QED', that is, processes involving fermions as well 
as bosons with QED-like interactions between them, and also the case where 
the bosonic field has a three-boson interaction: for these cases, almost-exact 
asymptotic results are easily found. The more complicated cases of QCD 
(with additional four-point boson vertices), scalar QED and 'all-out' theo- 
ries with arbitrary bosonic self-interactions, are discussed in section 4. 

2 Single-field theories 

The enumeration of tree diagrams in a given theory is simpler than that of 
general diagrams with loops, in the sense that for tree diagrams symmetry 
factors do not occur: indeed, at present the counting of higher-order dia- 
grams appears almost hopeless. Since diagrams can be counted by simply 
replacing all propagators and vertices by unity, the counting problem be- 
comes equivalent to solving the Schwinger- Dyson (SD) equation for a simple 
zero-dimensional Euclidean theory. To set the stage, let us consider such a 
theory, with a single self-interacting scalar field and a Euclidean action 

S{^) = \^'-V{ct>) , V{^) = ^^' + ^^' + --- (1) 

where V{(f)) collects all higher-point interactions according to whether they 
are present in the theory or not. Note that the minus sign and the factorials 
in front of the interaction terms ensure that all vertices are precisely unity. 
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Let us denote by a{n) the number of tree diagrams entering in the 1 — n 
amphtude, and by (f){z) its generating function: 



n>l 



^) = E rT«(^) • (2) 



This generating function, then, obeys the algebraic Schwinger-Dyson equa- 
tion 

<Piz) = z + V'i<Piz)) , (3) 

where double-counting is avoided by the factorials in front of the terms in the 
interaction potential: for instance, in pure ip^ theory the Schwinger-Dyson 
equation, translated back in terms of a(n), reads 

1 n\ 

a{n) = 6n,i + o E — j— 7^('^i)'^("'2)^«,"i+"2 • (4) 

^ ni,2>l ^1-^2! 

Obtaining the exact number of graphs for given n is simply a matter of 
algebraically iterating Eq.(H) up to the appropriate order in powers of z, 
which is a trivial task for any halfway decent computer algebra system. For 
the asymptotic result, however, we have to employ more. The asymptotic 
form of a{n) is of course given by the singularity structure of (f){z) as a 
function of z. Now, (f){z) cannot have poles for finite z if V{ip) is a finite 
polynomial, since the SD equation cannot then be satisfied. The singularities 
must therefore be branch points. Let us write z as a function of (p: 

z = F{<P)^<P-V'{<P) . (5) 

We now look for that value of (and z) for which the definition of (f){z) 
becomes ambiguous, i.e. where a branch cut starts. Such points 0o will be 
situated where d(f)/dz duverges, or 

dFU) 

= at = 00 • (6) 

To each of these roots corresponds a value zq = -F(0o), and that Zq which lies 
closest to the origin determines the leading asymptotic behaviour. We can 
then make an expansion around the appropriate value of 0o: 

z = F(0o) + V(0o)(0-0o)' + --- , (7) 
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and we can read off the approximate form of ^(z) in the neighbourhood of 
the singular point: 



1/2 



4>{z) - 00 - 1 - 



-2F((/.o) 



F{<f>o)J \ F"{<t>o) ■ 
The expansion of the square-root form, 

(2n-2)!a;" ^ 



n> 



\ n\{n - l)!22"-i 



E 



3/2^4^ 



(8) 



(9) 



where we indicate its asymptotic form, then tells us that the asymptotic form 
of a{n) is given by 



a{n) 



n 



3/2 



C 



1 



F{(t>o) 



D 



\ 27rF"(0o) 



(10) 



Two points are in order here. If it should happen that F" vanishes together 
with F', we might have to look for a cube-root branch point rather than a 
square-root one: wc have never yet encountered this case. Secondly, it is 
possible, as for instance in pure ip'^ theory, that there are several Zq values 
equally far from zero. In that case we have to add the asymptotic contribu- 
tions from every such point, and this is the mechanism by which in pure </7^ 
theory even values of n become forbidden. 

Some results arc collected in the the following table, for theories in which 
all interactions from if^ up to and including ip^^ are present. 



M 


C 


D 


3 


2.00000 


.282095 


4 


2.50804 


.191409 


5 


2.57845 


.178231 


6 


2.58755 


.175794 


7 


2.58859 


.175393 


8 


2.58868 


.175336 


9 


2.58869 


.175331 


oo 


2.58869 


.175329 



The last theory, with potential V{ip) = exp(93) — (p — 1, is interesting in that 
it establishes an upper bound on the number of tree graphs in any single-field 
theory. Equivalently, it gives the number of graphs arising from an effective 
action after tadpole renormalization. 
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3 QED-like theories 



We now turn to the more complicated case where fermionic fields are also 
present: we then have fermions, antifermions, and bosons. The simplest case 
is that of QED with a single fermion type. The action is then given by 

S{V,X,X) = '^^^ + XX-XX¥^ , (11) 

where x denotes the fermion and x the antifermion field. We now have 
three amplitudes, depending on the incoming line, and we have the following 
generating functions: 

(j){z,x,x) = — i — i — r(^i^ ^ noip,nix,n2x) , 
ip{z,x,x) = — i — i — r(^ix ^ noip,nix,n2x) , 

"0,1, 2>0 U 1 z 

- z"'''x^^x"'^ 

ip{z,x,x) = — i — i — r"'('^ ^ ^0<^,nix,n2x) , (12) 
not^>o nolniW 

and coupled SD equations: 

(j) = z + ipip , ip = X + , ip = X + (ptjj . (13) 

These can be readily expressed in alone: 

^ = , ^ = , (14) 

^ 1-0 ' ^ 1-0 ' ^ ^ 

and 

(^ = z+ (1^0)2 ' ^ = ^^ ■ (15) 

The combination ^ implies, of course, conservation of fermion number, and we 
see that it suffices to determine as a function of z and ^, except for processes 
without any external bosons. Again, combined expansion in powers of ^ and 
z is trivial by iterating Eg. (]T3|) . For the asymptotic behaviour we now have 
to study a two-variable problem. We do this by momentarily keeping z fixed, 
and considering the singularity of in terms of ^ : 

e = F,(0) = (1- 0)2(0 -z) , F;(0) = (1-0)(1 + 2z-30) . (16) 
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The singularity, parametrized by must be again of branch-point type, and 
is situated at 

1-1-2?' 4 
00 = ^ , F^{<Po) = ^{l - zf , F;(0o) = -2(1 - ^) , (17) 

so that the asymptotic behaviour for high powers of ^ is given by 

This immediately gives the behaviour with powers of z as well. Expanding 
the form (1 — 

A.(, A v.n,. J27/4r (3n-2 + fc)! 
gives for the number of graphs 

, (27/4)"(3n-2 + A;)!n!2 ^ ^ 

The goodness of this asymptotic result, when compared with the exact enu- 
meration, does not depend on k but only on n: the exact result is 1.3644356 
times the approximate one for n = 1, which ratio decreases to 1.0244771 for 
n = 10 and to 1.0120180 for n = 20. This /c-independence is related to the 
fact that we have here the exact Taylor expansion of the pole around z = 1: 
when we use Stirling's approximation for the factorials, a dependence on k 
is introduced. 

We may extend our discussion to the case of more fermions. If we intro- 
duce / fermion flavours, each flavour j will have its own generating functions 
ipj and ipj, with variables Xj and xj: but (j) will still be described by Eq.(|5]), 
with the sole redefinition 

/ 

^ = XX ^ ^ = ■ (21) 

The amplitude therefore becomes, upon multinomial expansion: 



a{(p kip, niXiXi, ^2X2X2, • • • , nfXfXf) 

(27/4)" (3n - 2 + ky.nlriiW ■■■rifl 
^ x/27^ (3n-2)!n3/2 



(22) 
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For theories in which the various fermion types have the same type of vertices, 
this is the general way in which one goes from the single-fermion to the multi- 
fermion case, and therefore we shall only consider the single-fermion case in 
the following. 

The next simplest case is that where we allow a ip^ self-interaction for the 
boson in addition to the V^XX vertex. The analog of Eq. ([T5|) is then 

so that 

= (1- 0)2(0 -l^^-z) (24) 

The equation F^{(j)) now has two roots, leading however to the same result: 
with = 1 — $ and z = 1/2 — (, the singularity condition reads = (, and 
the form of close to the singularity reads 

In turn, this gives the asymptotic form 



23n+fc-5/2^2n + A;-3/2)! 



n 



|2 



a(0 ^ A;0, nxx) - ,3/2^ (2n-3/2)!" " " ^^'^ 

Again, the goodness of the approximation is independent of k, and the ratio 
exact /approximate reads 1.253314137 for n = 1, 1.019251423 for n = 10, and 
1.009498692 for n = 20. 



4 Asymptotics by saddle-point methods 

It is tempting to extend the reasoning of the previous section to more com- 
plicated cases. This is, however, dangerous because of the following reason. 
In finding the singularity we have first to solve F^{(j)o) = to find 0o as a 
function of z, and then to determine the coefficient of z'' in 1 / Fz{4>o{z))"^ , 
where m is a large number. In principle, this is then again determined by 
the precise nature of the singularity of this form, that is, the behaviour of 
Fz{4>o{z)) close to a root. If zq is this root, we can generally write 

FziMz))=Aizo-zr{l+Bizo-z)+0{iz,-zy)) , (27) 
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for some positive p. Naively making the expansion around z — Zq gives then 



mB 



k 



+ 



ipm + k-l)\ ^ (pm + k-2)\ 

— mBzQ^-rr, ^ttt- + 



= E 



-) 



z \^ {pin + k - 1)! 



k\{pm — 1)! 



A^z^yzoJ k\{pm-l)\ 



Bzo 



k\{pm-2y. 

m{j)ni — 1) 
{pm + k — 1) 



+ 



(28) 



The correction term is not small for generic large m and k values, but only 
becomes small in the 'super-asymptotic' limit k ^ ^ 1, which is too 
asymptotic to interest us here. 



The most efficient way to determine the high-A; behaviour appears to be 
the following. Let us rewrite the relation between 0o and z: 



F'^{ct>) = ^ G{ct>) = z 



(29) 



For instance, in the 'QED+c/?^' case considered above, we have G(0) — 2(f) — 
0^ — 1/2. The form of the coefficient Cn,k of ^'^z''n~^/'^ is then given by a 
(counter-clockwise) Cauchy integral around z — 0: 



-FM) 



, G{ct>)^z 



(30) 



We can readily transform this into a loop integral for 0: 



1 

2m 



-H{4>) 

G{(pf+^H{(t)r\ 27riJ2(0) 
'g(</>)( 



Fg(^)(0) ,i/2(0)=F^(,)(0) 



(31) 



Note that, in the definition of i72, the substitution z = G{(f)) must be made 
after the differentiation. The only tricky point is to determine which of the 
various roots of G{4>) = should be chosen to integrate around. Note that 
this question can be answered unambiguously: if worst comes to worst, one 
can simply check the result against the exact enumeration for largish values 
of n and k. Having found the point around which to do the (J) integral, we 
then proceed to deform the integration contour into a pair of integrals parallel 
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to the imaginary axis. The upwards integral ('main' integral) is chosen to 
run over the saddle point ^, situated on the real axis and given by 



K'($) = , K{(l)) = -n log H{(j)) - k log 



(32) 



In many cases (j) can actually be given as a function of n and k in closed form. 
The downwards integral ('return' integral) is chosen to run over another 
saddle point. In the cases we have studied, we have always found that the 
values of H((f)) and -F(0) in the saddle point of the return integral are larger 
(in absolute value) than those for the saddle point of the main integral, so 
that the contribution of the return integral is exponentially suppressed with 
respect to the main integral. The result, therefore, is 



G'{4>) 



27rG'(0)fc+ii7(0)"\ i72(0)^"(0) 



(33) 



The final asymptotic estimate for the number of tree graphs with k + 1 
external bosons and 2n fermions/antifermions can therefore be written as 
follows: 



a{ip k(p,nxx) ~ 
C^ = l/H{^) , 
C2 = l/G(0) , 



D = 



-i7(0) 



27r \h2{^)K"{^) ■ 
The numbers Ci, C2 and D only depend on the ratio k/n. 



(34) 



As a first check, we redo the QED case with one fermion flavour. Here, 
G{(f)) = {3(f) — l)/2 so that we must integrate around = 1/3. There is only 
a single saddle point = (n + k)/{3n + k). Since > 1/3 and main integral 
is indeed upwards, and the return integral can be moved out to infinity. The 
result for Cn,k is nothing but the Stirling approximation of the 'exact' result. 

The next case is that of 'QED+(/9^'. The equation G(0) = has two roots, 
1 ± y^l/2, of which 1 — ^1/2 is on the appropriate Riemann sheet. This can 
be seen from the fact that the saddle point = {2n+k — \/2ti? + nk) / {2n+k) 
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is to the right of this root, and has a positive value for choosing the 

other possible saddle point (with a + sign before the square root) leads to a 
negative G"(0) and consequently a negative asymptotic estimate. Again, the 
saddle-point result boils down to the Stirling approximation of the 'exact' 
one. 

A more complicated case is that of QED with a pure (f^ interaction added. 
We find 

G'(</)) = -1(50^-20-6) , 

^2(0) = l(l-0)(502-20-6) . (35) 

Note that the roots of H2 are also roots of H or G': this is a general occur- 
rence. The main-integral saddle point is a root of the equation 

(5n + 3k)(f)^ - 3{n + k)(f)'^ - 6(3n + k)(f) + 6(n + /c) = . (36) 

By standard methods, we can find the three roots of this equation: 

(2r7r 1 



0r -- 



5n + 3k 



2(3171^ + SOnk + 7ey/^sml — - -u\+n + k 



f 29n^ + YBn^k + 63nk'^ + 17k^\ 

u — arcsm ^ — — — „, ^, ^ ,^ , (37) 

V (31n2 + 30nA; + 7A;2)3/2 j ' ^ > 

with r = 0, 1, 2. The choice r = interpolates smoothly between .325259493 
for A;/n = and 1 for n/ A; = 0, and this turns out to be the correct saddle 
point. The singularity structure of the integrand is as follows: there are 
poles of order A; + 1 at the three roots of G, -1.788306912, .3252594930, and 
2.063047419. There is a pole of order 3n — 1 at 1, and poles of order n at 
±\/2. With the standard convention that the square root branch cut lies 
along negative real values, there are cuts along the real axis from — -\/2 

to (1 - a/3I)/5 .9135528726 and from (1 + \/3T)/5 ~ 1.313552873 to 

^/2. The return integral can be taken to cross the real axis at -.9135528726, 
where H = 2.041476117 and G = -1.344004635. The saddle-point values 
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for the main integral are always smaller than these in absolute value, con- 
firming the above statement that the return integral may safely be neglected. 



The more relevant case of QCD is treated in the same manner. We have 
G(0) = -^0^-^0^ + 20-1 , 
G'(0) = l(4-50)(2 + 0) , 
H{cp) = 1(1-0)3(2-20-02) , 

^^2(0) = -l(l-0)(4-50)(2 + 0) , (38) 

and solves 

{5n + 3k)(j)^ + 3(3n + k)(f)'^ - 12(2n + k)(f) + 6{n + k) = ; (39) 
this root can be written as 

' ^ 2{13P + 49n^ + 50nA;)^/2 sin (^^v}j - 3n - k , 




, Uln^ + 225n^k + 123nk^ + 23k^\ . . 

u = arcsin 2 ; — -— — . (40) 

' (13A;2 + 49n2 + 50nA;)3/2 J ^ ' 

The singularity structure resembles that of the previous case: poles of order 
A; + 1 at the three values -3.343142188, .2853836802, and 1.257758508; a pole 
of order 3?2 — 1 at 1, and poles of order n at —1 — -\/3 ~ —2.732050808 and 
at — 1 + -\/3 ~ .7320508076; and branch cuts running from —1 — to -2 
and from — 1 + -\/3 to 4/5. The loop-integration contour is situated around 
.2853836802 and moves smoothly from this value upwards to .7320508076 
as k/n increases from to infinity. The saddle point for the return integral 
is at = —2, where H = 27/2 and G = —25/6, again always considerably 
bigger than H{(f)) and G{(f)). This allows for the complete determination of 
Ci^2 and D: the value of is given by Eq.(^), and 

50^ + 1203 + 202- 360 + 20 
^'^^ = ^ (2-20-02)2(1-0)2 

250^ + 6003 ^ 5402 _ 2040 + 156 

(503 + 902 - 240 + 6)2 ■ ^ ' 
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In the following table we present the non-universal quantities Ci, C2 and D 
for various ratios k/n. 





'^i 




n 
u 


-o.U 


o.i4o 


4i00. 


.uuuzyzy 


-2.5 


8.157 


1316. 


.0005196 


-2.0 


8.217 


417.5 


.0009172 


-1.5 


8.389 


133.6 


.001594 


-1.0 


8.961 


43.78 


.002640 


-0.5 


10.89 


15.40 


.003802 





18.00 


6.452 


.003882 


0.5 


48.97 


3.666 


.0002183 


1.0 


197.2 


2.833 


.0006782 


1.5 


836.8 


2.599 


.0001506 


2.0 


3141. 


2.535 


.00002913 


2.5 


10660. 


2.516 


.000005336 


3.0 


34530. 


2.511 


.0000009583 



As expected, the accuracy of the asymptotic approximation improves uni- 
formly if n and k grow with a fixed ratio. In the table we collect some 
results, where we have of course only those values for which both n and k are 
integers, and have iterated the exact generating function up to n + A; = 29. 
The accuracy is actually quite reasonable even for moderate values of n and 
k. 
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n 


k — n/2) 


k = n/2 


k = 2n/3 


k — n 


k = 3n/2 


k^2n 


k — 2)n 


1 












1.019 


1.003 


2 




.9904 




1.010 


1.010 


1.007 


.9995 


3 


.9709 




.9999 


1.006 




1.004 


.9993 


4 




.9927 




1.004 


1.004 


1.003 


.9994 


5 








1.003 




1.002 


.9995 


6 


.9837 


.9946 


.9992 


1.002 


1.003 


1.002 


.9995 


7 
8 




.9958 




1.002 
1.002 


1.002 


1.001 
1.001 


.9996 


9 


.9888 




.9993 


1.001 




1.001 




10 




.9966 




1.001 


1.001 






11 








1.001 








12 


.9915 


.9971 


.9994 


1.001 








1 








1 nni 

l.UUl 








14 




.9975 




1.001 








15 


.9932 




.9995 










16 




.9978 












17 
















18 


.9943 


.9980 












19 
















20 
















21 


.9951 















We may also consider a theory where tp'^ bosonic self-interactions occur 
for every p. We simply list the results: 



30 - -0exp(0) - - exp(0) , 
1 

3 -exp(0) - -0exp((/)) , 

l(_l + 0)3(_2 + exp(0)) , 

-(-l + 0)(-6 + 2exp(0) + (/)exp(0)) , 
exp(0)((n + + n - A;) - (6n + 2A;)0 + 2k , 

12 — 10 oxi)(o) + 3oxi)(2o) — lo(xp(p) + 2o^(xi)(o) 
(-2 + exp((/)))2(-l + 0)2 

^(60exp(0) + exp(20) + 6</)2exp(0) + 36 - 24exp((/))) 
(-60 + 0exp(0) + exp(0))2 



G(0) = 

G'(0) = 

H{ct>) = 

= 

K"{<t>) = 



13 



The Riemann sheet structure of the function G is of course much more com- 
phcated in this case, but fortunately the relevant saddle point is the simplest 
solution on the real axis, again interpolating smoothly between the appro- 
priate zeroes of H and G. In the table, we give the non-universal quantities 
as a function oik/n. 



logio(A;/n) 


Ci 


C2 


D 


-3.0 


8.150 


2000. 


.0002904 


-2.5 


8.170 


1250. 


.0005151 


-2.0 


8.224 


384.6 


.0009093 


-1.5 


8.396 


135.1 


.001580 


-1.0 


8.977 


44.05 


.002614 


-0.5 


10.91 


15.53 


.003753 





17.99 


6.494 


.003797 


0.5 


48.08 


3.731 


.002092 


1.0 


184.5 


2.907 


.0006325 


1.5 


729.9 


2.681 


.0001372 


2.0 


2564. 


2.618 


.00002613 


2.5 


8333. 


2.591 


.000004754 


3.0 


20000 


2.581 


.0000008515 



These numbers are qualitatively quite similar to that for QCD. The accuracy 
of the approximation appears to be almost identical to that of the QCD case. 

Another interesting case is that of scalar QED, where we have an addi- 
tional ^p^xX interaction term. The SD equations now become more compli- 
cated: 

V' = x + V'(0 + 072) , V' = ^ + V'(0 + 072) , = -z+(l + 0)V'V' , (43) 
so that 

Following the same steps as before, we find 

-2 + 902 + 403 + 60 



G'(0) = 



3(2 + 20 + 02) 
4(02 + 20 + 4)(l + 0)2 
3(2 + 20 + 02)2 



14 



(2 - 20 



12(2 + 20 + 02) ' 

(02 + 20 + 4)(-2 + 20 + 02) 



2 + 20 + 02 

the main-integral saddle point is a solution of 

{An + A;)0^ + 3(3n + A;)0^ + 6n0 - 2(n + /c) = 

and reads 

1 



(45) 



(46) 



= 



u 



in + k 
arcsin 



2(A;2 + 4nA; + 71^)^/2 sin 

n(3A;2 + 18n/c + 25n2) \ 



71 + U 



3n — k 



(A;2 + Ank + ^2)3/2 
Finally, we also need 

4n(8402 + 260^ + ( 



(47) 



^"^"^^ (-2 + 20 + 02)2(2 + 2 

8A;(36 + 32002 + 3710^ + 2050^ + 720^ 



+ 06 + 24 + 560 + 6403) 



+ 



k2\2 



+ 



+ 



1480 + 43203) 



(2 + 20 + 02)2(-2 + 902 + 403 + 60)2 



(48) 



The results for the non-universal constants are given in the following table. 



logio(A;/n) 


Ci 


C2 


D 


-3.0 


9.671 


4189. 


.1731E-5 


-2.5 


9.690 


1326. 


.5466E-5 


-2.0 


9.756 


420.2 


.1721E-4 


-1.5 


9.970 


133.7 


.5371E-4 


-1.0 


10.66 


43.18 


.1627E-3 


-0.5 


13.06 


14.55 


.4459E-3 





23.03 


5.510 


.8568E-3 


0.5 


85.11 


2.667 


.6268E-3 


1.0 


805.8 


1.776 


.1574E-3 


1.5 


.1526E5 


1.496 


.2247E-4 


2.0 


.4021E6 


1.407 


.2542E-5 


2.5 


.1197E8 


1.379 


.2647E-6 


3.0 


.3712E9 


1.370 


.2681E-7 
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The non-universal constants appear to vary much more rapidly as a function 
of k/n than in the case of QCD. The accuracy of the asymptotic estimate, 
however, is essentially the same. 

Conclusions 

We have demonstrated how the number of tree-level diagrams for several 
theories can be computed exactly for given numbers of external legs of various 
kinds. We have also described how asymptotic formulae for these numbers, 
valid in the limit of many legs of each kind, can be obtained. We have 
compared these results for several theories of interest, including QED, QCD 
and scalar QED. Comparison with the exact results up to fairly high order 
shows that the asymptotic estimates are accurate. 
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